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ABSTRACT 

The authors's point of view, based on QCD, on the nuclear quark structure 
is presented. Different models for explaining the EMC-effect are considered. 
It is also shown that cumulative production data are very useful for a bet- 
' ter understanding of the EMC-effect and give some evidence in favor of its 

■ multiquark nature. 

Discovery of the EMC-effect || has drawn attention of the world-wide community 
i ! of physicists to the problem of quark structure of nuclei, and to its irreducibility to the 
quark structure of constituent nucleons only 0. Stream of theoretical papers followed 
the discovery of EMC suggesting a whole spectrum of possibilities for understanding the 
phenomena 0]. However, many of the suggestions met with difficulties after a change 

■ of experimental data on F a /Fd in the region of small x Nowadays, when all 
suggestions seem to be made, one can try to analyze them on a general basis and to 

Q-i! estimate to what extent the nuclear quark structure is understood and what is still unclear. 

1. Connection of nucleon and nuclear quark structure. Probably 
^ . G.West first noticed |7j that QCD evolution equations result in a simple convolution 
J> ! relation of nonsinglet quark distribution functions (the valence quarks) of nucleus and 

\ nucleon [] 

H : A 

xF 3A » V A (x,Q 2 ) = jf daT% a (a)V N (^,Q 2 ) (1) 
where the function T A s (a) satisfies the baryon number sum rule 

daT A ws (a) = l (2) 

(all nuclear function here are normalized to A). Due to this, one can consider T A (a) as 
an effective "valence nucleon" distribution function over a fraction of momenta a in spite 
of the impossibility of expressing it through a one-nucleon wave function. 

1 This is an updated version of the talk jy at 13 International Seminar on Rclativistic Nuclear Physics 
and QCD (Dubna, Sept. 2-7, 1996) which includes new experimental data and some ideas appeared since 
that time. Partialy supported by the RFBR Grant 96-02-17631. 

2 Since the evolution equations do not depend on the kind of object 

V A (n, Q 2 )/V A (n, Q 2 ) = V N (n, Q 2 )/V N {n, Q 2 ) = j N (a s (Q 2 )) 

(the dot means derivative with respect to log Q 2 and n is the number of a moment) the first equality 
immediately gives V A (n, Q 2 ) = T^ s (n)Vjv(n, Q 2 ) which is equivalent to (pi). 



o 



I 



A similar relation can be written for the singlet channel as well || [| which mixes the 
singlet quark, T,(x, q 2 ) = J2 q x[q(x, Q 2 ) + q(x, Q 2 )} ~ F 2 , and gluon, G(x, Q 2 ), distribution 
functions 



Z A (x,Q 2 )= [ daT^(a)E N (-,Q 
Jx \a 

G A (x,Q 2 ) = ( A daT s A (a)G N (-,Q 2 ) 

J x \ Oi J 



(3) 
(4) 



where, in general, Tf ^ T A and Tf satisfies the energy-momentum sum rule 



da aT^(a) 



M A 
AM, 



N 



(5) 



Really, diagonalizing the system of two linear evolution equations for the moments 
S J 4(n, Q 2 ) and Ga(^, Q 2 ), one can obtain the relation for two eigenfunctions f ± (n, Q 2 ) = 
Q 2 ) + C^(as(Q 2 ))G(n,Q 2 ) (C^ are some diagonalizing coefficients depending on 
anomalous dimension matrix elements): 



fi(n,Q 2 ) 



T±(n)f±(n,Q 2 ) 



(6) 



The quark (and gluon) distribution function is expressed through the limit of quark 
(or gluon) propagator (P|g(0)g(£)|P) when £ — > 0, regularized with the help of an ul- 
traviolet cutoff parameter Q 2 . It must satisfy Bethe-Salpeter equations (Fig. 1) with 
inhomogeneous terms. The second equation for the gluon propagator however become 
homogeneous, at least in the so-called leading logarithm approximation, i.e. disregarding 
a^(logQ 2 ) n ~ 1 corrections. In this approximation also, it became an algebraic equation for 
the moments of structure functions with coefficients which are independent of a target. 
This means that the ratio E(n)/G(n) is also independent on the target and as a result of 
(^) T\ = = Tjf which leads to the relations (|3|,[|) . 

Physically this approximation corresponds to a rather widely accepted scheme in which 
the gluon distribution is totally due to QCD evolution process starting from the valence 
quark distribution at low Q 2 , i.e. when there is no "primordial" gluon distribution (see 
e.g. the work JIU[). 
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Fig. 1. Bethe-Salpeter equation for singlet quark and gluon structure functions. 



An immediate consequence of Eqs. (Q^) is the equality of average momenta fractions 
of gluons and quarks in the nucleus and nucleon 



< XG >A _ < Xq >A 

< Xg >N < X q >n 



(7) 



This relation is in good agreement with BCDMS [[J data which are the most precise 
nowadays, e.g. < x q >^ 2 / < x q >d 2 — 1 = (0.7 ± 1.7 ± 1.0)%. The difference of average 
momenta for other nuclei is also zero within the error bars (see Table 2 in Ref. fT0|j). 

The relation (0) clearly contradicts the very popular rescaling hypothesis [|ll]] in ex- 
planation of the EMC-effect. In fact, the passage from nucleon to nucleus in these models 
is equivalent to the growth of Q 2 for which, according to QCD, < xq > increases and 
< x q > decrease^. 

In conclusion of this section let us stress once more that QCD evolution equations 
just as relation (||) are results of the leading twist approximation. So, the relations ([[]) 
and (0^) do not include the nuclear screening which is, at least formally, a high-twist 



effect ||12| , [L3|. Some experimental observation of a significant Q -dependence of Fs n /Fc 
in the region x < 0.02 was known recently [TTJ . 

2. The EMC— effect. Let us see now what the EMC-effect means in the frame of 
our approach. Let us assume that the functions Ta determine, at least approximately, an 
effective distribution of nucleons in nucleus and therefore they are mostly concentrated in 
the region of a = 1 (i.e. in the region of zero internal momentum of a nucleon). Expanding 
F/v(^) in (|I]) and (|3|) around a = 1, it is easy to obtain for not very large x 



R = jr 



t< T > + < (1 



F' 1 

a)T > x^ + - < 
' F N 2 



a f T > x ( x ^L + 2^] + 
V F N F N I 



where < • ■ • > means integration over interval < a < A. If one accept that F^ ~ 
(1 — x) k and k ~ 3, then x-dependence of the second and the third terms are the factors 
— k[x/(l— x)] and k[x/(l — x)] ■ [(k — l)x/(l — x) — 2] respectively. In the region of x ~ 0.5 
the second term is close to zero and to obtain the depletion of R from unity in this region 
discovered by EMC one should have 



< T 



-1 = A, > 



or 



da 



T b Aa)-T 



NS, 



a) 



for the ratio R 2 of the structure functions F 2 ~ £ and 



1- < aT? s >=S A >0 



or 



da a 



a 



T 



NS, 



a 



A A >0 



5 A > 



(9) 



(10) 



for the ratio R3 of the structure functions of xF^. 

In addition, in the region x ~ 0.5 the sea quarks are practically absent: therefore one 
can expect here R3 ~ R 2 and 



5 A ^ A A 



(more exactly 



The relations @ and (0), mean that the number of "effective nucleons" in a nucleus 
has to be larger than A, and the valence nucleons have to carry only a part of the total 



'Another criticism of the hypothesis from QCD point of view can be found in Ref. 
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nucleus momentum. In other words, the EMC-effect is the result of a repumping of a part 
of momentum from valence quarks to sea quarks in the nucleus in comparison with free 
nucleon. 

Notice, that the shock produced by the discovery of EMC was due to the prejudice 
that a nucleus is made of A nucleons and so the condition A^ = has to be imposed 
on the distribution T s , which unavoidably results in R^{x ~ 0.5) = 1, independent of 
the form of T s . In this sense, the difference between T s and T NS (necessary to explain 
the EMC-effect) leads to the irreducibility of the nuclear quark structure to the quark 
structure of free nucleons. 

In spite of generality, this approach allows one to draw a number of interesting con- 
clusions: 

i) It immediately follows from ([J) that the ratio f] 

R 2 (x~0)= daT^ia) = 1 + A A > 1 (12) 
Jo 

The most accurate measurement of BCDMS [Ej shows a small ~ 5% but definite excess 
of the ratio over 1 in the region of small x, i.e. the same value as the loss of momenta 
of the valence nucleons 5 A . This means a small number of particles of the non-nucleon 
component which have to be heavy enough to supply the 5% repumping of the momentum 
(e.g. p-mesons, iViV-pairs or pions far off the mass shell). 

ii) In addition to the internucleon sea there is a small, ~ A^, but hard enough "col- 
lective sea" of quark-ant iquark pairs in nuclei. 

Using ([[]) and (|) it is easy to obtain for the sea 

o a (x,q 2 ) = j: a -v a 

= f A daT% s {a)0 N (-,Q 2 ) + [* da [Tf (a) - T% s (a)} S w (-,Q 2 ) (13) 

where the first term comes from the internucleon sea, which rapidly decreases with in- 
creasing x, and the second term comes from a "collective sea", 0' A , which is hard since 
its center of gravity is 

ao,= <n- T r> = a~a 1 (14) 

For pions on the mass shell this number is m^jm^ m 1/7. That is the reason why the 
repumping of the momentum into the pions [|T^] gives no satisfactory description of the 
data in the region of small x (too many pions are needed to supply the 5% repumping). 

in) The place of intersection i^O^o) = 1 does not depend on the sort of nucleus and is 
at xq ~ 0.3. Really, if there are no screening and light particles in nuclei, T A (a) has to 
be smooth enough in the region of small x. Using then the first two terms of (H) for R 2 
it is easy to find 

x _ 1 ( _ J X ° da aTf (a) \ ~ 



1 — xo 3\ Jo° daT A (a) J 

The ratio of integrals in the r.h.s. is in the interval [0, xq] and thus 0.25 < xq < 1/3. This 
feature of the ratio R 2 has been well confirmed experimentally |16| with xq = 0.278±0.010. 



'Recall that the screening phenomena are disregarded here 
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Now, what about different models proposed? They are, in fact, different suggestions 
of the repumping mechanisms. Not all of them seem satisfactory from our viewpoint. We 



have mentioned the rescaling models [[II] where part of the repumping comes into gluon 
component. However, the main drawback of these models is the softness of the gluon 
and the sea component in nucleon. This leads to an extra big value of ^2(0) after the 
5% repumping. (Although the authors deny the applicability of their model to the region 
of small x.) As it was noticed, models with repumping of momenta into the mass shell 
pions |L5l have the same disadvantage. 



Other models can be divided into three big categories: 



i) Models with repumping of the momentum either into massive meson component |u]] 
(p, u, off the mass-shell pions) or into nucleon-antinucleon pairs 0. A component like 
that is probably related to the core of nuclear force at small distances. However, it is 
hard to believe that the nucleon can conserve at such small distances its individual quark 
structure without converting it into multiquark states. 

ii) Repumping inside each nucleon |18| , i.e. change of its quark structure due to the 



influence of the internuclear field. Transition of part of nucleons into A-isobar [19| also 
belongs to this class. We do not see, however, how it is possible to obtain the hard sea 
here. 



Hi) Repumping inside a multiquark fluctuation |20|. By this we mean not only a 



bound state but also a state of two or more nucleons with interaction of their quarks, 



as proposed in fl2fl , or with an exchange quark interaction in the final state considered 
in That kind of interactions is inevitable in any theory with a composite nucleon. 
However, the calculation of the quark structure of the states like that seems as difficult as 
the calculation of the quark structure of nucleus. Recently some progress in this direction 
has been achieved |22] . 



It is necessary to stress the important difference between a multiquark state and few- 
nucleon correlation (FNC) ||23|| . The loss of momenta of the valence quarks for the latter 
are the same as averaged over the nucleus, Afnc = Aa, due to a change of structure of 
each nucleon. For the multiquark, however, it has to be much larger 

A 6q > A A (15) 

e.g. if there is no repumping inside the nucleons, then A^ = paAq q , where pa is a 
probability of multiquark states. In fact, the relation (^) can be considered as a definition 
of the multiquark state in distinction with FNC. A statistical realization of the hard 



antiquark sea is also known (see e.g. Kondratyuk paper [2(J). 

It seems that structure function measurements alone cannot distinguish between these 
models. So, new sources of information are necessary. One of them is deep inelastic scat- 
tering with measurement of hadrons in a final state. Production of p- and A-resonances 
and also i^~-mesons and antiprotons which carry the information about the collective sea 
is especially interesting for evident reasons. 

3. Cumulative particles production. Another source of information is the 
cumulative particle production. Especially, the production of fT~-mesons and antiprotons 
off nuclei in the region x > 1, because of the peculiarity of the nuclear quark structure 
mentioned before. 

A question arises however: to what extent is the cumulative production cross section 
determined by the nuclear structure functions Fa(x)7 Until now there are no quite reliable 
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data for nuclear deep inelastic scattering in the region x > 1, though there are some 
indications of similarity of the cumulative meson spectra and structure function F 2 (x) in 
this region |25| . 

There exist two points of view on the physics of cumulative production ||: (a) "hot 
models" , in which massive clusters in nuclei (which are necessary to produce a cumulative 
particle) are formed by an incoming hadron either by a sort of compression of the nuclear 
matter and heavy fireball formation or by multiple rescattering ||26|| ; (b) "cold models", in 
which formations of that sort already exist in nuclei because of Blokhintsev's fluctuations of 
density |27| either in a form of multiquark states or in a form of a few nucleon correlation, 



resulting in the high momentum tail of Fermi motion. This reflects in the structure 
functions of the nucleus. A common property of these models is the independence of the 
nuclear parton fragmentation of the nucleus type. This allows us to write down the cross 
section of the process in the form 

= PA->h(x,y,Pr) = I — F A (a)f h (-,y,p T ) (16) 
A d 6 p Jx a \a / 

where x = —u/s, y = —t/s and the function fh does not depend on A, i.e. it is the same 
for a nucleus and for a free nucleon. Combining (|IIJ) with ([TLB!), it is easy to obtain a 
natural expression 

pA^h(x,y,PT) = J daN A (a)p N ^ h (^,y,pT^+J daN A (a)p^ h (^,y,p T ^ (17) 
where we use the notation 

N A = \ (T| + Tr) and N A = \ (t| - T^) . (18) 



The first expression in (|18|) can be considered approximately, due to smallness of the 
EMC-effect, as a distribution of nucleons over fractions of the momentum. For cumulative 
and stripping protons it is necessary to add to (|T7| ) a term proportional to N A (x) which 
takes into account dissociation of the nucleus. Moreover, just this term gives a major 
contribution when px ~ [28]. Parametrizing the form of the spectrum of stripped and 
cumulative protons with px ~ (with normalization < N A >= 1 + A A /2, < aN A >= 
1 —5 A /2) and using the experimental cross section for pN-*n, one obtains the cross section 
of cumulative pion production without any new parameter. (The second term in ( P~T| ) 
naturally gives a small correction.) This program for deuterium (to minimize possible 
secondary nuclear effects) has been made in work [PJ and shows a good agreement with 
experiment. Also, the ratio K + /tt + agrees with experiment. This agreement confirms the 
independence of fragmentation of the kind of a nucleus (at least, for light nuclei), which is 
the base of (|TT|,|i8|) and means also that the valence mesons carry the same information on 



the nuclear quark structure as the cumulative protons pO ]. However, the main peculiarity 
of the nuclear quark structure is hidden here. 

Interpretation of N A in (|17| , |T8"|) depends on the mechanism of repumping and, due to 
the second term in (|TTP, dominates for "sea particles" K~,p in the region x > 1. They are 
just sensitive to the peculiarity of the nuclear quark structure. For the ratio of K + /K~ 
yields in the region we have 

_K^ J x A daN A (a) PN _ tK+ ^) , . 

A K~ ~ f x A daN A (a)p N ^ K+ ^) ' 
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where the approximation p^^ K _ ~ Pn->k+ is used. It is known experimentally |3T], 
p3| that the ratio for various nuclei (Be, Al, Cu, Ta, Pb) is constant in x, within 
experimental accuracy, in the region 1 < x < 2.5 (Figs. 2) [] is in sharp contradiction 
with its behaviour in region x < 1. Therefore, the functions Na and A"^ in this region 
should only differ by a factor. One can assume that this difference reflects the difference 
in normalization condition for this functions 

< aN A >~< Na >= A A /2 and < aN A >~< N A >~ 1 (20) 

So for the models of type i) and ii) one would expect from fll9|) pa 2/ A a- Using the 
parametrization [| of the SLAC [|4j and BCDMS § data for the EMC-effect one finds 
r Be = 86 (A Be = 0,023), r Al = 55 (A Ai = 0.036), r Ta = 36 (A Ta = 0.056) and r Pb = 35 
(Apt, = 0.058), which is significantly higher than the experimental ratio (see Figs. 2), 
especially for light nuclei. 
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Fig.2. The relative yields of cumulative K+ and K~ from works |l) (Ei ab = 8.9 GeV, 9i ab = 168° 
(10 GeV, 119°) and @ (40 GeV, 159°). 



For the repumping inside multiquark states, which 
has to determine the cumulative cross sections in this 
region of x, the repumping A Qq ~ A a/ Pa has to be larger 
(due to a small pa) and ~ 2/A 6(J has to be lower. E.g. 




the experimental ratio r Be ~ 9 corresponds to A§ q ~ 
0.22 and p^e — 0.10. So, the low ratio of can be 
considered as an indication of the multiquark mechanism 
of both cumulative phenomena and the EMC-effect. 

Now, let us turn to the cumulative antiprotons. Natu- 
rally, they are sensitive to the A^A^-pair repumping mech- 
anism 0. The ratio of p/p yields is determined by 
an expression of the type of (0) and is of the order 
2/Aa — 10 2 . The experiment [[JIJ gives for this ratio 
p/p pa 10 8 , which definitely rejects the above mechanism. On the other hand, with no 
packing of the collective sea into AW-pairs, the cumulative p can arise in fragmentation 
of q (just as K~). Then the ratio p/K~ has to be ~ 0.1 (suppression by an order of 



Fig.3. 
works 



2.6 X 



The K /p yields ratio from 
H {10 GeV, 119°). 



(xp 



Notice that the variable x here is defined by a minimal mass (in the nucleon mass units) of M| = 
a + Pb — Pc) 2 for the process B + A —* C + X, where pa is a 4-momentum of a nuclei per a nucleon. 
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magnitude is due to fragmentation q — > p and some growth due to a smaller transverse 
momentum of p at the same value of x) which is not far from the experimental observa- 



tion |35| K~ jp 5 for Be (see Fig. 

More accurate calculations of K + /K~ andp/ K~ ratios based on the dual string model 
with the multi-quark fluctons were done in paper [|36|| . Notice, however, that the model 
underestimates the K + / K~-iaXio for the proton target in the region x > 0.6. 

It is necessary to stress also that secondary nuclear effects can be significant for the 
intermediate and heavy nuclei we have considered. Indications to these effects come, for 
example, from an enhanced A-dependence of cumulative proton and K + productions and 
from an enhanced by 4-5 times depletion from unity of pa-^-kI Pd-*-k i n the region x ~ 0.6 
compared to the deep inelastic scattering. (Notice also that the ratio of cumulative 
cross sections He/ D shows even an anti-EMC effect in this region!) For these reasons it 
would be very desirable to have accurate data on the kaon and antiproton production off 
deuterium. 

The conclusive headlines are: 

i) The cause of the EMC-effect is the repumping of the valence quark momentum to 
a collective sea of quark-ant iquark pairs. 

ii) Small excess of the A/D ratio in the x ~ region points to a hardness of a collective 
sea or to a massive non-nucleon component in nuclei. 

Hi) Many popular models are in trouble due to i) and ii). 

iv) The ratio of K + /K~ cumulative cross sections supports the multiquark mechanism 
of the EMC-effect and of the cumulative process. 

v) Accurate data on the kaon and antiproton production off deuterium are highly 
desirable. 
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